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ABSTRACT 

We model the X-rays reprocessed by an accretion disk in a fiducial low-mass X-ray binary system 
with a neutron star primary. An atmosphere, or the intermediate region between the optically thick disk 
and a Compton-temperature corona, is photoionized by the neutron star continuum. X-ray lines from 
the recombination of electrons with ions dominate the atmosphere emission and should be observable 
with the Chandra and XMM-Newton high-resolution spectrometers. The self-consistent disk geometry 
agrees well with optical observations of these systems, with the atmosphere shielding the companion 
from the neutron star. At a critical depth range, the disk gas has one thermally unstable and two stable 
solutions. A clear difference between the model spectra exists between evaporating and condensing disk 
atmospheres. This difference should be observable in high-inclination X-ray binaries, or whenever the 
central continuum is blocked by absorbing material and the extended disk emission is not. 

Subject headings: accretion, accretion disks — atomic processes — instabilities — line: formation — 
X-rays: binaries 



1. INTRODUCTION 

The X-rays reprocessed on accretion disks to the opti- 
cal, UV and X-ray bands in both X-ray binaries and AGN 
could be one of our best probes of accretion disk structure 
and the environment around compact objects. It is con- 
venient to study an illuminated disk around neutron star 
low-mass X-ray binaries (LMXB) since the central contin- 
uum dominates the energetics of the disk atmosphere and 
accretion disk corona (ADC), and the system geometry 
is often determined. Measurements of the optical mag- 
nitudes and lightcurve amplitudes of LMXB yielded esti- 
mates of the angle subtended by the disk of ~ 12° and of 
the average fraction of X-rays reprocessed by the disk into 
the optical of < 0.1 (de Jong, van Paradijs, & Augusteijn 
1996). Models that seem to correctly describe the accre- 
tion disk interior have been unable to produce disk geome- 
tries consistent with the optical data (Dubus et al. 1999). 
Analytical models assuming an isothermal disk show bet- 
ter agreement with the data, yet produce illuminated disks 
that are significantly thinner than expected (Vrtilek et al. 
1990). 

A pure continuum incident on the gaseous disk will pho- 
toionize and heat the gas. Radiative heating sustains a 
~ 10 7 K ADC and a smaller ~ 10 5 -10 6 K atmosphere emit- 
ting X-ray lines on top of the optically thick disk. X-ray 
spectra of LMXBs during dips show strong Comptonized 
emission plus an excess flux at 0.65 keV consistent with 
the same covering fraction (Church et al. 1998), positing 
that ADC and disks may coexist. The apparent size of 
the ADC and the disk may not be equal (Church 2000), 
since the emission mechanisms have different radial depen- 
dence. Nayakshin, Kazanas & Kallman (2000) modeled 
the reprocessed continuum and fluorescence emission from 
the inner disk of AGN with an ADC, and Ko & Kallman 



(1994) modeled X-ray spectra of disk annuli in an LMXB 
with ADC and normally incident radiation. 

In this paper, we calculate the structure of an accretion 
disk illuminated by the grazing radiation of a neutron star 
in an LMXB with a Compton-temperature ADC, building 
on a model by Raymond (1993). This model, which we will 
detail in a future paper, includes a more accurate calcula- 
tion of the vertical and radial disk structure than previous 
models, accounts for thermal instabilities, and incorpo- 
rates recent atomic model calculations that predict X-ray 
spectra at a resolution superior to the Chandra and XMM- 
Newton spectrometers. We obtain the first self-consistent 
flared disk geometry that matches the disk thickness ex- 
pected from optical data, and we predict that broad X-ray 
lines will be observable and will provide, among other di- 
agnostics, evidence for photo-evaporation, which is central 
to the origin of winds (Woods et al. 1996). 

2. MODEL ATMOSPHERE 

We consider an LMXB with a M* =1.4 M© primary 
radiating an Eddington luminosity bremsstrahlung contin- 
uum with T = 8 keV. We use a set of fiducial system pa- 
rameters for a bright LMXB, so application to a particular 
source will require using the observed X-ray continuum to 
improve accuracy. The maximum radius of the centrally- 
illuminated disk is 10 11 cm, so the orbital period ~ 1 day. 
The minimum radius is 10 8 ' 5 cm, below which the omitted 
effect of radiation pressure will dominate. 

We obtain the vertical disk atmosphere structure of each 
annulus by integrating the hydrostatic balance and 1-D ra- 
diation transfer equations for a slab geometry: 
dP_ _ GM^pz 

dz r 3 
dF v _ k v F v 

dz 
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while satisfying local thermal equilibrium: 

A{P,p,F„) = (4) 
where P is the pressure, p is the density, F v is the incident 
radiation field, Fj} is the reprocessed radiation propagat- 
ing down towards the disk midplane, z the vertical dis- 
tance from the midplane, G the gravitational constant, 9 
the grazing angle of the radiation on the disk, r the radius, 
and k„ is the local absorption coefficient. The difference 
between heating and cooling A includes Compton scatter- 
ing, brcmsstrahlung cooling, photoionization heating, col- 
lisional line cooling, and recombination line cooling from 
H, He, C, N, O, Ne, Mg, Si, S, Ar, Ca, and Fe ions. We 
assume cosmic abundances (Allen 1973). The code from 
Raymond (1993) computes the net heating and ionization 
equilibrium, Compton scattering, and line scattering us- 
ing escape probabilities. A new disk structure calculation 
simultaneously integrates Eqs. f-3 by the Runge-Kutta 
method using an adaptive stepsize control routine with er- 
ror estimation, and Eq. 4 is solved by a globally convergent 
Newton's method (Press 1994). At the ADC height z cori 
the equilibrium T is close to the Compton temperature 
T compton-, from which we begin to integrate downward un- 
til T < Tphotir)- The optically thick part of the disk, with 
temperature T p hot, is assumed to be vertically isothermal 
(Vrtilek et al. 1990). To get Tphot, we assume that the 
viscous energy and illumination energy is all locally radi- 
ated as blackbody radiation. That is, for z p hot *C r and 
R* < r: 

4 3GM,M (l-y)L x sinfl(r) 

^Phot - 8?rr3 + 7^ IP) 

where i?* is the neutron star radius, M, L x the accretion 
rate and luminosity, a the Stephan-Boltzmann constant, 
and r\ — 0.9 is the X-ray 'albedo' derived from optical 
observations (de Jong, van Paradijs, & Augusteijn 1996), 
such that 1 — 77 is the fraction of X-rays absorbed at the 
photosphere. The height at which the integration ends 
is defined as the photosphere height z p hot- Thus, the as- 
sumption is that for z < z p hot viscous dissipation domi- 
nates heating. 

The boundary conditions at the ADC are set to 
P(zcor) PcorkT com pt on / piTip, F x (z cor *) — Lj x j^KT , and 
F^(z cor ) = 0, where F x = J F v dv, k is the Boltzmann 
constant, and p is the average particle mass in units of the 
proton mass m p . The boundary conditions at z p hot for F v 
and F^ are set free, and we use the shooting method (Press 
1994) with shooting parameter p cor adjusted until we sat- 
isfy P(zphot) = PphotkTphot/ prrip at the photosphere. Note 
Pphot is the viscosity (a) dependent density calculated for 
an X-ray illuminated Shakura & Sunyaev (1973) disk. 

We calculate 9(r) iteratively. To get T p hot self- 
consistently from Equation 5, we need 9(r) = 
arctan(dz atm jdr) — arctan(z atm /r) + arctan(i?*/r). We 
neglect the term oc R*/r, which is valid for r > 10 8,5 
cm. After an initial guess for z a t m (r), we define it as the 
height where F x is attenuated by 1/e. We stop after 9{r) 
and T p hot converge to < 10 %. Note F x — > for z < z p hot- 
Since z atm is not physically determined, but it is bound by 
Zatm > Zphot , we ran a second model calculating 9(r) from 
Zphot, which we will use to estimate systematic errors due 
to the 1-D radiation transfer calculation. 



3. THERMAL INSTABILITY OF THE PHOTOIONIZED GAS 

The disk atmosphere exhibits a thermally unstable re- 
gion, which has direct consequences for the disk structure 
and X-ray spectrum. The origin of the thermal instabil- 
ity is well understood, and its behaviour depends on metal 
abundances and the local radiation spectrum (Hess, Kahn, 
& Paerels 1997). 

Within a range of pressure ionization parameters (S), 
thermal equilibrium is achieved by three distinct temper- 
atures, only two of which are stable to perturbations in 
T (Fig. 1). 5 ee 4irP rad /P gas , where P rad = F x /4wc is 
the radiation pressure and P gas the gas pressure (Krolik, 
McKee & Tarter 1981). The thermally stable branches 
have positive slope (Field 1965). Previous spectral studies 
of illuminated accretion disks in LMXB had not explic- 
itly selected the stable solutions (Ko & Kallman 1994), a 
choice which affects X-ray production. 

The instability implies a large temperature gradient as 
the gas is forced to move between stable branches, requir- 
ing the formation of a transition region whose size may 
be determined by electron heat conduction, convection, or 
turbulence, depending on which dominates the heat trans- 
fer. For simplicity, we neglect emission from the tran- 
sition region. Upon calculation of the Field length Xp, 
the lengthscale below which conduction dominates thermal 
equilibrium (Begelman & McKee 1990), we estimate that 
conduction forms a transition layer <~ 10~ 2 times thinner 
than the size of the X-ray emitting zones. Nevertheless, 
X-ray line emission from the neglected conduction region 
may not be negligible in all cases (Li et al. 2001). 

Conduction determines the dynamic behaviour of the 
gas. The dynamics will allow us to attach a physical in- 
terpretation to the chosen solutions from Fig. 1 . Roughly, 
a static conduction solution can only be found if the tran- 
sition region splits the instability S-range in half (Fig. 1). 
This solution (Zeldovich & Pikelner 1969) takes the low-T 
stable branch at S < S stat , and the high-T stable branch 
at S > S stat , separated by the transition layer at E sta t- 
A transition layer located away from S stat will dynami- 
cally approach S stat by a conduction driven mass flow. A 
transition layer at E evap > E sta t is moving towards lower 
S, and the gas is evaporating from the low-T branch to 
the high-T branch. A transition layer at £ co „<j < S stat is 
moving towards higher S, with the gas condensing to the 
low-T branch (Zeldovich & Pikelner 1969; Li et al. 2001). 

We compute the disk structure for both condensing and 
evaporating solutions. We assume a steady state, condens- 
ing or evaporating mass flow through the transition layer 
at Econd or E evap , respectively. The static conduction so- 
lution is an intermediate case of the latter extreme cases. 
We take a single- valued T(S), since a two-phase solution 
would be buoyantly unstable, making the denser (colder) 
gas sink. The evaporating disk corresponds to the low- 
T branch, while the condensing disk corresponds to the 
high-T branch (Fig. 1). This introduces spectral differ- 
ences (section 6). 

We do not know from first principles whether the disk 
atmosphere is evaporating, condensing, or static. We 
estimate the speed of the conduction mass flow to be 
Vcond = 2kT '/iPgasXp, by using the characteristic conduc- 
tion time at the Field length, where n is the Spitzer (1962) 
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conductivity (for a detailed discussion, see McKee & Begcl- 
man 1990). We obtain a conduction mass flow speed v cone i 
which is 1-2 xl0~ 2 times the local sound speed. Thus, the 
phase dynamics will depend on the subsonic (v > v con d) 
flow patterns in the disk atmosphere, and these flows will 
also determine the evaporation or condensation rates. 

The only physical mechanism known to transport the 
necessary angular momentum for disk accretion involves 
a magneto-rotational instability (MRI) which drives tur- 
bulent flow in the disk (Balbus & Hawley 1998). Heat 
transfer due to this turbulent flow could quench the ther- 
mal instability and affect the disk structure. The MRI 
also favors our assumption of vertical isothermality in the 
optically thick disk. 

4. SPECTRAL MODELING 

With the disk structure p{r,z), T(r,z) and ion abun- 
dances fz,i{f,z), we model the X-ray line emission from 
the disk atmosphere by using HULLAC (Hebrew Univer- 
sity/Lawrence Livermore Atomic Code, Klapisch et al. 
(1977)). The code calculates the atomic structure and 
transition rates of radiative recombination (RR) and ra- 
diative recombination continuum (RRC) for the H-likc and 
He-like ions of C, N, O, Ne, Mg, Si, S, Ar, Ca, and Fe XVII- 
Fe XXVI. The emissivities are calculated as described in 
Sako et al. (1999). 

The spectrum for each annulus is added to obtain the 
disk spectrum. Each annulus consists of a grid of zones 
in the vertical z direction, and T, p and fz.i+i for each 
zone are used to calculate the RR and RRC emissivities. 
The radiation is propagated outwards at inclination an- 
gle i, including the continuum opacity of all zones above. 
Compton scattering is not included in the calculated spec- 
trum because it is negligible in the vertical direction in the 
region where the line emission is formed. The spectrum 
is Dopplcr broadened by the projected local Keplerian ve- 
locity, assuming azimuthal symmetry. 

5. DISK STRUCTURE RESULTS 

The self-consistent disk is thicker than would be ex- 
pected from the local pressure scale height Hp alone. We 
find Hp < Zphot < z a t m - We fit the modeled photo- 
sphere and atmosphere height, z p hot and z atm , with Cr n 
(r in cm), with fit parameters C and n. We find C p hot — 
2.4 ±0.4-1.9 x 10~ 3 cm 1 -"'", n phot = 1.14 ± 0.01+0.06, 
C atm = 1.0±0.1+0.2xl0- 3 cm 1 - n ««™,n o t m = 1.21±0.01. 
Thus, Zphot ~ 3-4 Hp and z atm ~ 7-8+2 Hp. Estimated 
systematics are shown, if significant. Vrtilck et al. (1990) 
estimated n atm = 9/7 = 1.29, but in spite of the steeper 
radial dependence, the Vrtilck disk is thinner, and it equals 
Hp for r > 10 10 cm. For comparison to the optically de- 
rived disk thickness, z a tm is the de facto disk boundary, 
since a fraction 1/e of the central X-rays have been ab- 
sorbed there. Thus, previous theoretical studies severely 
underestimated the disk thickness. 

The atmosphere's (z > z p hot) maximum photoelectric 
opacity is always t<1, although r/ sm8(r) > 1. Illumi- 
nation heating dominates at r > 10 10 cm (Eq. 5), where 
only F x (z p hot) / Fx{zcor) ~ 0.12 of the incident photons 
reach z phot directly, while F^(z phot )/(F x (z cor ) sin 6) ~ 
0.35 reach z p hot after reprocessing in the atmosphere. 
Thus, the atmospheric albedo is ~ 0.5. To reproduce the 



disk (atmosphere and photosphere) albedo observed by de 
Jong, van Paradijs, & Augusteijn (1996) (n <~ 0.9), the 
photosphere albedo can be < 0.8, which is a range closer 
to physical expectations. There is no significant change in 
the atmosphere structure by varying a from 1 to 0.1. 

6. SPECTRAL RESULTS 

The LMXRB photon flux (photon cm" 2 s" 1 kcV" 1 ) is 
modeled by: 

ptot = e -*„N* H p* + e -v„N$° k pdisk ^ 

where F* is the neutron star continuum, F. ^ lsk is the RR 
line and RRC modeled flux, N* H and N§ sk are the neu- 
tral hydrogen absorption column densities and a v are the 
Morrison & McCammon (1983) absorption cross sections. 
The system is assumed to be 10 kpc away. 

If we let N* H = Nfj lsk in Eq. 6 with i = 0, the lines are 
swamped by the continuum. Thus, low inclination neu- 
tron star LMXBs are unlikely to have detectable X-ray 
lines from the disk. Instead, we consider = 5 x 10 22 
cm~ 2 and Nfj sk = 10 21 cm~ 2 with i = 75°, where an ob- 
scuring medium absorbs half of the continuum flux and is 
compact enough to leave the disk almost unobscured. 

With a partially obscured central continuum, disk evap- 
oration has an observable spectral signature. We simulate 
50 ks XMM-Newton RGS 1 and Chandra MEG +1 obser- 
vations (Fig. 2,3). Some bright lines are shown on Table 
1. The evaporating and condensing disks have contrasting 
O VII/O VIII line ratios. The evaporating disk contains 
gas at T <~ 7-10 x 10 4 K, unlike the condensing disk. The 
H-like ion line intensities are higher for the condensing disk 
since it has more gas at T ~ 10 6 K. The spectral differences 
stem from the distinct differential emission measure distri- 
butions d(EM)/dlogE and from the O VII recombination 
rate app oc T~ 7 , where 7 = 0.7-0.8, EM — J nifi e dV is 
the emission measure, and n i7 n e are the ion and electron 
densities. 

The RRC and He-like line triplet can be used as tem- 
perature and density diagnostics, respectively. The RRC 
width oc T (Liedahl & Pacrels 1996). The forbidden line / 
of O VII at 22.097 A is absent due to collisional depopula- 
tion at high density, since p > 10 13 cm~ 3 . The O VII inter- 
combination i to resonance r line ratio is g = (/ + i)/r ~ 
i/r > 4, indicating a purely photoionized plasma (Porquct 
& Dubau 2000). The O VII RRC broadening is Doppler- 
dominated, resembling the RR lines, while the O VIII RRC 
has a component with FWHM ~ 2 A that is produced at 
T - 10 6 K. 

7. CONCLUSIONS 

Using a model for a disk illuminated by a neutron star, 
we have produced a global spectrum of X-ray line emis- 
sion from a self-consistent accretion disk atmosphere. The 
flaring of the atmosphere matches optical observations (de 
Jong, van Paradijs, & Augusteijn 1996). For fiducial pa- 
rameters of a LMXB with a partially obscured central 
source with Eddington luminosity, observable spectral sig- 
natures would identify an evaporating disk atmosphere. 

The O VII/O VIII line ratios are a good tracer of evap- 
oration from the disk. The line-rich spectrum has line 
profiles resolvable by the XMM and Chandra spectrom- 
eters. The lines have broad wings and a double-peaked 
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core with a velocity from the outermost radius. Only a 
single line may be resolved since the emission is well dis- 
tributed throughout 10 8,5 cm < r < 10 11 cm. He-like 
triplet diagnostics show signatures of photoionization and 
high density (n e > 10 13 cm~ 3 ), and H-like RRC are very 
broad, with a superposed narrow component in evaporat- 
ing disks. 

The distinct spectral signature of evaporation is due to 
the presence of a thermal instability in the photoionized 
gas of the disk. This instability cannot be quenched by 
spectral variations, abundance variations, or optical depth 
effects (Hess, Kahn, & Paerels 1997). Conduction heat- 
ing should not affect this result, though turbulent heat 
transfer or coronal heating might quench the instability if 
present. In the outer disk, the incident radiative power 
is > 10 times the viscous power, so coronal-like heating 



should be small. 

The X-ray lines from the disk are weak and difficult to 
detect for a low-inclination LMXB. Only if the inclina- 
tion is high and a cloudy environment exists do we expect 
them to be detectable. At press time, an XMM-Newton 
observation of LMXB dipper EXO 0748-67 has detected 
broad lines from photoionized gas with n e > 7 x 10 12 cm~ 3 
(Cottam et al. 2001), which we interpret as the disk at- 
mosphere emission, since it confirms the key predictions of 
our model. 

We are grateful to S. Kahn and I. Joseph for fruitful 
discussions. This work was performed under the auspices 
of the U.S. Department of Energy by the University of 
California Lawrence Livermore National Laboratory un- 
der contract No. W-7405-Eng-48. 
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log E = log (F x /cP gas ) 

Fig. 1. — Ionization parameter H vs. temperature assuming an 8 keV brcmsstrahlung continuum. The s-curve corresponds to the locus of 
solutions in thermal balance and ionization equilibrium. We model the extreme cases, the evaporating and condensing disks, which correspond 
to the lower and upper stable branches of the curve, respectively. The inclusion of conduction in the energy equation modifies the thermal 
equilibrium locus, creating solutions with one small transition region. Schematics of the static transition region at E. B tat and two dynamic 
transition regions at H conc j and H et)ap are shown. For logH > 1.2, T increases up to T comp ton ~ 10 7 K (not shown). 
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Fig. 2. — Spectra for evaporating and condensing accretion disks. Simulated 50 ks observation with XMM-Newton RGS 1. The continuum 
emission from the inner (r < 10 8,5 cm) disk is not included, and the obscuration of the ionizing continuum from the compact object is 
assumed. 
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Fig. 3. — Condensing disk spectrum (as in Fig. 2). Simulated 50 ks observation with Chandra MEG 1. 

Table 1 

Line fluxes for two disk models. 



Line(s) Evaporating 51 Condensing 51 



O VII Hea triplet 


2.1 


0.99 


O VIII Lya 


3.4 


4.7 


O VII RRC 


1.3 


0.61 


O VIII Ly/3 


0.94 


1.5 


Fe XIX L (0.815 keV) 


0.45 


0.91 


Ne X Lya 


1.4 


2.2 


Fe XXVI Ha (1.28 keV) 


0.74 


0.86 


Mg XII Lya 


1.2 


2.0 


Si XIII Hea triplet 


1.2 


2.1 


Si XIV Lya 


4.0 


6.0 


Fe XXV Lya 


15 


16 


Fe XXVI Lya 


12 


12 



a Line fluxes in units of 10 -12 erg cm -2 s — 1 , typically over 
an interval E/AE ~ 100. Estimated systematic normaliza- 
tion error limit of ~ 50 %, due to the 1-D calculation. 



